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Abstract 

We construct the most general physically admissible positive-definite inner product on the 
space of Proca fields. Up to a trivial scaling this defines a five-parameter family of Lorentz invari- 
ant inner products that we use to construct a genuine Hilbert space for the quantum mechanics 
of Proca fields. If we identify the generator of time-translations with the Hamiltonian, we ob- 
tain a unitary quantum system that describes first-quantized Proca fields and does not involve 
the conventional restriction to the positive- frequency fields. We provide a rather comprehensive 
analysis of this system. In particular, we examine the conserved current density responsible for 
the conservation of the probabilities, explore the global gauge symmetry underlying the conser- 
vation of the probabilities, obtain a probability current density, construct position, momentum, 
helicity, spin, and angular momentum operators, and determine the localized Proca fields. We 
also compute the generalized parity (V), generalized time-reversal (T), and generalized charge 
or chirality (C) operators for this system and offer a physical interpretation for its VT-, C-, and 
CT^T-symmetries . 

1 Introduction 

Quantum Mechanics (QM) and Special Relativity (SR) constitute the basis for the main body of 
modern physics developed during the first half of the twentieth century. The conceptual marriage 
of QM and SR is Relativistic Quantum Mechanics (RQM) in which the one-particle quantum wave 
functions are identified with the solutions of an appropriate field equation. The solution space of these 
equations provide the representation (Hilbert) spaces for the irreducible (projective) representations 
of the Poincare group that in turn define the elementary particles [U |21 [3] . One can construct the 
Fock space associated with such a Hilbert space and define the quantum field operators as operators 
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acting in this space. In this sense, RQM provides a basis for Quantum Field Theory (QFT). Although 
a fully relativistic treatment of elementary particles requires the methods of QFT, there are processes 
and energy scales where one can safely employ RQM. This is for example the case in the relativistic 
treatment of quantum systems displaying quantum superposition and quantum entanglement [1]. 

Since the early days of RQM the formulation of a consistent (first-quantized) quantum theory 
with a genuine probabilistic interpretation for bosonic fields has encountered severe difficulties. It 
was in an attempt at such a formulation for the massive scalar (Klein-Gordon) fields that Dirac 
introduced the electron's wave equation and discovered the antimatter. His proposal that led to the 
development of second quantized field theories was also mainly motivated by this problem. Recently, 
we used the methods of Pseudo-Hermitian QM [5] to give a complete formulation of a genuine 
quantum mechanical treatment of Klein-Gordon (KG) fields [6] . The purpose of the present paper is 
to employ a similar approach to devise a consistent quantum theory of free first-quantized massive 
vector fields. 

Research on the theory of massive vector fields started with the work of Proca [7J whose original 
aim was to obtain a description of the four states of electron-positron system using a four-vector. 
Proca attempted to pursue the approach of Pauli and Weisskopf [H] who had quantized the KG fields 
and interpreted the KG conserved current with that of the electric charge rather than the probability. 
Although Proca's results did not serve its original purpose, its mathematical formalism could be used 
to treat massive vector fields. 

In its manifestly covariant form, Proca's equation for a free massive vector field A^^ reads [7] 

d^.F^"' - M^A"" = 0, (1) 

where 

F^" := d^'A" - (2) 

^ ._ mc -g ^j^g inverse of the Compton's wave length, m G M"^ is the mass, d^d^ := rj^^d^du-, and 
{vn = M-^ :=diag(-l, 1,1,1). 

If we apply to ([1]) and use F^^, = —F^^, we find that A^ satisfies the Lorentz condition: 

d.A'' = 0. (3) 

Making use of this relation in ([1]) we find 

{d^d'' -M^)A^' = 0. (4) 

This means that the four components of the vector field A^ satisfy the KG equation. Note that the 
Proca equation does not have a gauge symmetry. Therefore ([3]) is not a gauge choice but rather a 
constraint that is to be imposed on A'^. In fact, it is not difficult to show that (jl]) together with 
([3]) are equivalent to the Proca equation ([1]). Similarly to the case of KG fields, the presence of 
the second time-derivative in (jlj) is responsible for the difficulties associated with devising a sound 
probabilistic interpretation for the Proca fields. 
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An important motivation for the study of Proca fields is their close relationship with electro- 
magnetic fields. For example, the study of Proca fields may shed light on the interesting problem 
of the consequences of a nonzero photon rest mass which has been the subject of intensive research 
over the past several decades (See the review article [9].) It might also be possible to construct a 
first-quantized quantum theory of a photon by taking the zero-mass limit of that of a Proca field. 
This can in particular led to a resolution of the important issue of the construction of an appropriate 
position operator and localized states for a photon \10\ . 

The literature on the Proca fields is quite extensive. There are a number of publications that deal 
with the issue of the consistency of Proca's theory and the difficulties associated with interacting 
massive vector bosons (See [HI [12] and references therein.) It turns out that there are various ways 
of formulating a relativistic wave equation describing the dynamical states of a massive vector boson. 
The famous ones employ the equations of Proca, Duffin-Kemmer-Petiau (DKP) [13] , and Weinberg- 
Shay-Good jH]. These are equivalent in the absence of interactions but lead to different predictions 
upon the inclusion of interactions. 

Among the works that are directly relevant to the subject of the present article are those of 
Taketani and Sakata [15j, Case [16], and Foldy [17J. Using the analogy with the Maxwell equations, 
Taketani and Sakata |15J reduced the ten-component DKP wave function to a six-component wave 
function satisfying a Schrodinger equation. They made use of this representation to study the 
interaction of the field with an electromagnetic field. Employing the six-component Taketani- Sakata 
(TS) representation. Case |16j cast the Proca equation into a form in which the positive- and negative- 
energy states were separately described by three-component wave functions. The latter is the so- 
called Foldy's canonical representation [17| which Case used to study the non-relativistic limit of 
the Proca equation and obtained the position and spin operators acting on the TS wave functions. 
The approaches of Case [T6j and Foldy [IT] involve the use of an indefinite inner product |18j on the 
space of six-component TS wave functions. This seems to be the main reason why these authors 
did not suggest any reasonable solution for the problem of the probabilistic interpretation of their 
theories. The problem of the construction of an appropriate position operator and the corresponding 
localized states for Proca fields has also been considered by various authors, e.g., [191 EOl EH [22] . 
But a universally accepted solution has not been given. Some other more recent articles on Proca 
fields are [231 [21]. 

In |25[ [26l [27l [H [28], we have employed the results obtained within the framework of Pseudo- 
Hermitian Quantum Mechanics (PHQM) to formulate a consistent quantum mechanics of KG fields. 
Here we pursue a similar approach to treat the Proca fields. The first step in this direction has been 
taken by Jakubsky and Smejkal [29] who constructed a one-parameter family of admissible inner 
products on the space of Proca fields. In what follows we will offer a more systematic and general 
treatment of this problem. In particular we give a complete characterization of Lorentz-invariant 
positive-definite inner products that render the generator of the time-translations and the helicity 
operator self-adjoint. We further construct a position operator and the corresponding localized states 
for the Proca fields and examine a variety of related problems. 

In the remainder of this section we briefly review PHQM and give the notations and conventions 
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that we use throughout the paper. 

PHQM [HI [30] has developed in an attempt to give a mathematically consistent formulation 
for the PT-symmetric Quantum Mechanics [3T]. In PHQM, a quantum system is described by 
a diagonalizable Hamiltonian operator H that acts in an auxiliary Hilbert space Ti' and has a 
real spectrum. One can show (under some general conditions) that these conditions imply that H 
is Hermitian with respect to a positive-definite inner product that is generally different from the 
defining inner product (■, ■) of [32j. This motivates the following definition of a pseudo-Hermitian 
operator [33]: H is a pseudo-Hermitian operator if there is a linear, invertible, Hermitian (metric) 
operator rj : Ti' ^ Ti.' satisfying = rjHrj~^ . This condition is equivalent to the requirement that 
H be self-adjoint with respect to the (possibly indefinite) inner product {{■,-))ri '■= "'?■)• It turns 
out that for a given pseudo-Hermitian operator H, the metric operator r) is not unique. If we make 
a particular choice for rj, we say that H is rj-pseudo-Hermitian. 

A proper subset of pseudo-Hermitian operators is the set of quasi- Hermitian operators |34j . 
A quasi-Hermitian operator is a pseudo-Hermitian operator that is ?7_|_-pseudo-Hermitian for a 
positive-definite metric operator 77+, i.e., 



This means that H is Hermitian (self-adjoint) with respect to the positive-definite inner product 



In particular if is a diagonalizable operator with a real spectrum and can be mapped to a Hermitian 
operator via a similarity transformation [32] . 

A quasi-Hermitian Hamiltonian operator defines a quantum system with a consistent probabilistic 
interpretation provided that one constructs the physical Hilbert space Ti of the system using the inner 
product ((% ■))^^ and identifies the observables of the theory with the self-adjoint operators acting in 



Although PHQM employs quasi-Hermitian Hamiltonian operators, the study of the pseudo- 
Hermitian Hamiltonians has proved to be essential in understanding the role of anti-linear symmetries 
such as PT . As shown in Ref. [35j, given a diagonalizable pseudo-Hermitian Hamiltonian if, one 
may introduce generalized parity (P), time- reversal (T) and charge or chirality (C) operators and es- 
tablish the C-, PT-, and CPT-symmetries of H that would respectively generalize C, PT, and CPT 
symmetries [36] of PT-symmetric quantum mechanics^] We recall that H is P-pseudo-Hermitian 
and T-anti-pseudo-Hermitian, V and T are respectively linear and anti-linear invertible operators, 
and C is a linear involution (C^ = 1) satisfying C = rj^^V, [H^lfl 

Next, we give our notations and conventions. 

^Ti. is the Cauchy completion of the inner product space obtained by endowing the span of the eigenvector of H 
(in Ti.') with the inner product ((•, •))r;+- 
2 See also [37]. 

simple consequence of the latter relation is that the CPT-inner product of f36^ is a particular example of the 
inner products ([6]) of [32 l [38 l l39] . 



(5) 




(6) 
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Throughout this paper a = (a°, a) stands for a four-vector a^, Greek indices take on the values 
0, 1, 2, 3, Latin indices take on 1,2, 3, and Sijk denotes the totally antisymmetric Levi-Civita symbol 
with £i23 = 1. We employ Einstein's summation convention over repeated indices and use (Tq and Aq 
to denote the 2x2 and 3x3 identity matrices, respectively. Recall that ctq together with the Pauli 
matrices 

^Ol\ f -i\ /^lO 

10' ^^^U ' '^^^ -1 



(7i = 



(7) 



form a basis of the vector space of 2 x 2 complex matrices, and that {ctq, cr^} is a maximal set of 



commuting matrices. We denote their common eigenvectors by e+ := 
the labels correspond to the eigenvalues of a^. 
Similarly, Aq together with 
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(8) 



form a basis for the vector space of 3 x 3 complex matrices, and {Aq, A3, Ag} is a set of maximally 



commuting matrices whose common eigenvectors we denote by e+i :— 
I \ 









: 













and 



eo := 




1 ) 



. Again the labels and ±1 are the eigenvalues of A3. 



Using the bases {(jj} and {Aj} we can construct the basis {Sm} for the vector space of 6 x 6 
complex matrices, where for all m e {0, 1, 2, • • • , 35} 



A., if 



m 



i + 4j. 



(9) 



The operators Eq, S3, E12, S15, E32, and E35 form a maximal set of commuting operators with 
common eigenvectors 

6^,5 := Ce (g) 65, (10) 

where e e {+, — } and s e {—1, 0, +1}. 

The organization of the paper is as follows. In Section 2, we present a six-component formulation 
of the Proca equation, establish its relation with pseudo-Hermiticity, and construct a positive-definite 
metric operator and the corresponding inner product. In Section 3, we compute the generalized 
parity 7^, time-reversal T, and charge or chirality C operators and elaborate on the C-, T'T-, and 
CT'T-symmetries of the Hamiltonian. In Section 4, we derive the expression for the most general 
physically admissible positive-definite inner product on the solution space of the Proca equation and 
demonstrate the unitary-equivalence of the representation obtained by the corresponding Hilbert 
space and the generator of time-translations as the Hamiltonian with the Foldy representation. In 
Section 5, we obtain and explore the properties of a conserved current density that supports the 
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conservation of the probabilities. In Section 6, we introduce a position basis and the associated 
position wave functions for the Proca fields. In Section 7, we discuss the position, spin, helicity, 
momentum and angular momentum operators, and construct the relativistic localized states for the 
Proca fields. In Section 8, we compute the probability current density for the spatial localization 
of a Proca field. In Section 9, we study the gauge symmetry associated with the conservation of 
probabilities. Finally in Section 10, we present a summary of our main results and discuss the 
differences between our approach of finding the most general admissible inner product and that of 
Ref. m. 



2 Pseudo-Hermiticity and Proca Fields 

2.1 Covariant Dynamical Field Equation and Helicity States 

To begin our investigation, we briefly formulate the covariant dynamical theory of Proca fields and 
review the physical polarization and helicity states. 

Consider a Proca field such that for all x° G M and n E {0, 1, 2, 3}, J^^ (fx x)p < oo. 

Then we can express (jl]) as a dynamical equation in the Hilbert space L^(R^) © L^(R^) © L^(]R^) © 
L^(M^), namely 

A{x^) + DA{x'')=0, (11) 

where a dot denotes a x°-derivative, for all /x G {0, 1, 2, 3} and G M the functions A^{x^) : M'^ ^ C 
are defined by 

A^(x°)(x) := A^(x°,x), Vx G M^ 
and D : L^(R^) L^(R^) is the Hermitian operator: 

[D/](x) := [m' - VV(x), V/ G L'{R'). (12) 

Clearly, for all /i G {0,1,2,3} and G M, A'^(x°) belongs to L^(R^), and D is a positive-definite 
operator with eigenvalues cj^ := k'^ + M'^ and the corresponding eigenvectors 0k(x) := {2n)^'^^'^e'^^'^ = 
(x|k), where k G R^, k := |k|, and (-I-) denotes the inner product of L^(R^). 

It is important to observe that a solution of ffTTl) is a Proca field provided that for all x° G R the 
functions A(x°) fulfill the Lorentz condition that we can express as the constraint: 

£[A(a;°)] := i°(a;°) + ■ A(a;°) = 0, (13) 

where for all / G L^(R^), (.^/)(x) := — iV/(x), and " ■ " is the usual dot product. 

The Lorentz condition (JT3i) implies that there are only three independent components of A; using 
this condition we can express A^{x^) and A^{x^) in terms of A(a;*^) and A(x'^). This means that the 
initial data for the Proca equation ffTTj) is given by (^A{xq) , A^Xq)^ for some initial value a;Q G R of 
Xo- We shall however employ a manifestly covariant approach and treat all components of the field 
yl(x°) on an equal footing. As a result we define the complex vector space (V) of solutions of the 
Proca equation (fTT!) as 

V := {A := (A°, A) | Vx° G R, [d^ + D]A^(x°) = 0, i:[A(x°)] = } . (14) 
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It is not difficult to show that the Lorentz condition f|T3l) may be imposed in the form of a constraint 
on the initial data for the dynamical equation flTT]) . nameljo 

= = 0. (15) 

This in turn implies 

V = |a := (A°, A) I Va;° G M, [dl + D]A^(x°) = 0, and 3x° G M, £[^(0;°)] = i:[i(a;°)] = } . 

(16) 

We can express the solution of (ITT!) in terms of (A'^(xg), A^(xg)) according to 



A'^(x°) = cos[(a;° - xl)D^'^]A^'{xl) + sin[(x° - xl)D^'^]D-^'^ A^'^xl). (17) 

Note that here and throughout this paper we use the spectral resolution of D to define its powers, 
:= ^3 d^k{e + ^2)«|k)(k| for all a G M. 
Next, consider the plane wave solutions of (j4]): 

^^(k,a;x) =iV,,k<(k,a)e-^^->'^Vk(x), (18) 

where N^^i^ are normalization constants (Lorentz scalars), and a(;'(k, a) denotes a set of normalized 
four- dimensional polarization vectors fulfilling the Lorentz condition, i.e., 

a,^(k,a)<(k,a') =5..', fc^<(k,a) = 0. (19) 

These relations show that the Proca field f|T8|) has only three physical polarization states (cr G 
{1, 2, 3}). In a fixed reference frame in which the plane wave has momentum k, we choose to work 
with a pair of transverse polarization vectors: [21 SI] 

a,(k, 1) = (0, a,(k, 1)), a,(k, 2) = (0, a,(k, 2)), (20) 

whose space-like components, ae(k, a) with a = 1,2, are normalized vectors perpendicular to k. As 
a longitudinal polarization vector we choose [21 SI] 

a,(k,3)= (-,-^^), k^ = eu;^, (21) 

which together with ([20 
pleteness relation [3], SI] 



form an orthonormal set. These polarization vectors also fulfil the com- 



1 

^<(k,a)<(k,a) =77^^^ + — A;^r. (22) 

(T=l 



"The in £[A(a;g)] stands for -DA{x°a). 

^The longitudinal polarization vector used in Ref. [55] is not normalized. Note also that, unlike the authors of 
Ref. [29j . we do not fix the transverse polarization vectors a^(\i, cr), cr = 1,2. Although we fix a reference frame in which 
their time-like components are zero, their space-like components are still arbitrary. Instead of fixing these vectors, 
we only fix the longitudinal polarization vector ([2T|) and make use of equations and ([22]) to do the necessary 
calculations throughout this paper. 
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It turns out that the use of the circular polarization vectors, 

u,,o(k) :=a,(k,3), M,,±i(k) := i= [a,(k, 1) ± ^ a,(k, 2)] , (23) 

simphfies many of the calculations. This is because their spatial components are eigenvectors of the 
helicity operator: 

where (i^i, i^a, -^3), ^ := -^/l^l, and 

S,:=Xr, 82:= -Xn, := X2, (25) 

are the angular momentum operators in the spin-one representation of the rotation group [42j. They 
satisfy 

SiSjSk + S^SjSi = SijSk + 6kjSi, [Si, Sj] = ieijkS]^. (26) 

A simple consequence of these relations is the identity f)'^ = f) which we shall make use of in the 
sequel. 

We can also construct the basic (plane wave) solutions of (jl]) that have a definite helicity. These 
are given by 

yl^_,(k; x°, x) = iV,,,(k) <,(k) e-^^^''^Vk(x), e = ±,H = -l, 0, 1, (27) 
where N^^hi)^ are normalization constants. 

2.2 Six- Component Formulation and Pseudo-Hermiticity 

In analogy with electromagnetism, we express the antisymmetric tensor F^'^ in terms of the E and 
B fields US]: 

pij e'^'^Bm, F^' =: E\ (28) 
Then the Proca equation ([1]) takes the form 

B = VxA, A = -E- V^o, 

= -M-^V -E, E = M^A + V xB. ^ ' 

Next, we let 7i := {Y : M'^ ^ ^ Y, < 00} denote the Hilbert space of vector fields, where for 
all Y, Z : ^ C^, ^ Y, Z^:= /^g ci^x Y(x)* • Z(x). If we eliminate B and we can express ^ 
as a set of dynamical equations in the Hilbert space 7i, namely 

A(x°) = -E(a;°) + a^-2DE(x°), E(x°) = DA(a;°) + DA(x°), (30) 

where D : 7i —>■ 7i is the Hermitian operator: 

[DY](x):= [p-S)2-i^2]Y](x) = V(V-Y(x)), (31) 
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Y G TtI is arbitrary, and we have used the identities 

[i^'YKx) = -V2Y(x), [{K ■ S)Y] (x) = V X Y(x), (32) 

which also imply 

^■[{^- S)Y] = 0, ^-[i^- SYY] = 0. (33) 
Next we express (l30l) as the Schrodinger equation (also known as the Taketani-Sakata equation 

mm) 

in^(xO) = ff^(x°), (34) 
where for all a;" G M the state vector \E' and the Hamiltonian H are defined by 

7 G M — {0} is an arbitrary constant having the dimension of length, and Hi and H2 are the following 
Hermitian operators that act in 7i. 

Hi := 7(D + 2)) + -f-^M-^M^ - ®) = 7(^^ + ■ S)^) + -f'^M-\D - ■ S)^) , (37) 
H2 := 7(^ + 2)) - 7"'>f"'(>f' - = 7(a^' + ■ S)') - 7"'a^"'(^ - ■ Sf) . (38) 

Note that one can invert the first equation in fl5Ul) (or take a x^-derivative of the second equation in 
f l30|) and use E = -DE) to obtain 

E(a;°) = -D-^ [m^ + (i^ ■ S)^] A(a;°). (39) 

In view of 0351) and 0391) . the six-component state vector \1/ is completely determined in terms of A 
and A. This means that the space of the state vectors is isomorphic (as a vector space) to the 
space of the initial conditions (A(a;Q), A(a;o)). Because of the linearity of the Proca equation, the 
latter is also isomorphic to the vector space of Proca fieldsj^ 

The six-component vectors \l/(a;'^) belong to C'^^Ti. If we endow the latter with the inner product 
(■, ■) defined by 

2 

i=l 



^Similarly to the two-component representation of the KG fields [25], the choice of the six-component state vector 
([35]) is not unique. Its general form is 'i^g(x'>) = 5(2^°)^' where gix*^) G GL(6,C). ([30ll is equivalent to the Schrodinger 
equation ilVi/ g(^xO)[x^) = i?g(a;0)^'g(2,o-)(a;'^) where Hgf^^o) ■— g{x'^)H g{x'^)^^ -t- ihg{x^)g{x^)~^ . The arbitrariness in the 
choice of g{x^) is related to a GL{6, C) gauge symmetry of the six-component formulation of the Proca fields. We will 
take g{x'^) to be the identity matrix. This is a partial gauge- fixing, because we do not fix 7. Changing 7 corresponds 
to a gauge transformation associated with a Gi(l,R) subgroup of Gi(6,C), namely 'i>^ ^l*^/ = 5(7, 7')^'^ where 

(7(7, 7') = ^ ^ "z + 'y ^ , ^ (8) Aq. This gauge symmetry has no physical significance, and as we shall see our final 

results will not depend on 7. 
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/ CM ~ 

for all ,^ = ( ^2 j ' ^ ~ ( ^2 1 ^ C^®7i, and denote the corresponding Hilbert space, namely 7i©7i, 

by 7i', we can view ^(x*^) as elements of Ti' and identify H with a linear operator acting in H.' . One 
can easily check that H : Ti' ^ Ti' is not a Hermitian operator, but it satisfies H"^ = SaifSg^, 
i.e., it is Ss-pseudo-Hermitian. This implies that H is self-adjoint with respect to the inner product 
((■) ■))s3 := ('5 S3 ■) on C?) H. This in turn induces the following inner product on the space f lTll) of 
the Proca fields. 



^ A, ( A' + V > < ( A + V , A' ^ 



(41) 



where g G is a constant. Because of Ss-pseudo-Hermiticity of if, this inner product, which is 
sometimes called the Proca inner product, is invariant under the time-evolution generated by if, but 
it is obviously indefinite. 

The dynamical invariance of (14T1) is associated with a conserved current density, namely [l3l [2] 



(x) := ig [A,(x)*F'^^(a;) - F^^(x)*A,(x)] . (42) 



This is the spin-one analog of the KG current density. We can express the Proca inner product in 
a manifestly covariant form in terms of the current density (H2l) . But because this inner product 
is indefinite, it cannot be used to make the solution space of the Proca equation into a genuine 
inner product space. Again in analogy to the case of KG fields, one may pursue the approach of 
indefinite-metric quantum theories [H] and restrict to the subspace of positive-energy solutions (e.g. 
see [B]). But this scheme has the same difficulties as the corresponding treatment of the KG fieldslH 
Next, we will use the positive-definiteness of D to show that if is a diagonalizable operator 
with a real spectrum. This suggests that it is ?7+-pseudo-Hermitian for a positive-definite metric 
operator 77+ [32j. Therefore, H is Hermitian with respect to the positive-definite inner product: 
((■5 ■))r;+ := {■,V+')- The construction of 77+ requires the solution of the eigenvalue problem for H and 
iit. 

The eigenvalue problem for the Hamiltonian fl36|l may be easily solved. H has a symmetry 
generated by 

A := do ® f), (43) 

which is the helicity operator in the six-component representation of the Proca fields^] The simulta- 
neous eigenvectors of H and A are given by 

*..(k):^if : " Tlk! r -kit! ) 

2 V ('^fe -erfc)u,,;,(k) +7r;. kM^,/^(k) / 
ii^,,;,(k) = f;,(k)^,,,(k), A^,,^(k) = h^,M, (45) 



^For example, although the inner product (j4ip restricted to the positive-energy solutions is positive definite, (x) 
that should correspond to the probability density is not generally positive-definite even for positive-energy fields. 
^Hereafter we will omit (X" wherever there is no risk of confusion. 
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where e G {-1, 1}, h G {-1, 0, +1}, k G M^ A; := |k| 



cuk := VP + A?, Tk ■■= a/o^, E^i^) := ehuj^, 

Ue,hO<) are circular polarization vectors ( l23l) . and 0^ are eigenvectors of D corresponding to the 
eigenvalues cuk- The eigenvectors "^^^hO^) together with 

form a complete biorthonormal system for the Hilbert space. This means that 

(^,,,(k),<l>,,,,(k')) = 5.,e'5/^,h'5'(k-k'), [ ^^'k|^,,,(k))(<l>,,,(k)| = So, (47) 

where (■, ■) stands for the inner product of Ti.', and for ^, C ^ 'H', \0{C\ is the operator defined by 
\0{(\x '■= (C;X)^; for all X ^ T^'- Similarly, one can check that indeed $e,/i(k) are simultaneous 
eigenvectors of W and A"!" = A with the same eigenvalues, W^^ hiX) = -E'e(k)$e,/i(k), A^^ h(k) = 
h $e,h(k), and that H and A have the following spectral resolutions 

H = YY I ^^'kE,(k)|*,,,(k))(<l>,,,(k)|, ^ = YY I ^^'k/i|^,,^(k))($,,,(k)|. 

Another remarkable property of the biorthonormal system {\l'e /i(k), $e /^(k)} is that 

(^,,;,(k),^,,,,(k')) = (^e',/.'(k'),^.,/^(k)) = ee'(<l>,,,(k),$,,,,,(k')). (48) 

Using the properties of the polarization vectors, i.e., Eqs. (IT^ and (122]) . we can compute the 
positive-definite metric operator, : 7i' — > 7i', associated with the biorthonormal system 
{^,,/,(k),$,,;,(k)}. The result is 



e=± h=0,±l 

The inverse of ?7+ has the form 



^+ = E E / rf'k|<i.,,,(k))(<i.,,,(k) 



H2 Hi 



(49) 



= E E Lrf^k|vl.,,(k))(vl.,,(k)| = ^ f ) . (50) 



e=± /i=0,±l 



Now we are in a position to compute ((■,■)),,+ . For all ^, C ^ "^'^ we let ^*,C\^±'C± G be 
defined by 

^^j:=e, (^^)^=^' ^±-^'±^'' C±:=C^±C^. (51) 
Then in view of (ET]), (ISH]), and ^ 
1 



2 



(52) 
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If we view Ti' as a complex vector space and endow it with the inner product fl52|) . we obtain a new 
inner product space whose Cauchy completion yields a Hilbert space which we denote by /C. 
Next, for all x'' G M we define U^o : V ^ 7i' according to 



t/,oA:=l./^vI/(a;°), VA G V, (53) 

where n G is a fixed but arbitrary constant. We can use U^o to endow the complex vector space 
V of Proca fields with the positive-definite inner product 

{{A A')) := {{UM,UM'h, = -^mx'),^'{x'))),,. (54) 

Because r/+ does not depend on the inner product ((■, ■))^^ is invariant under the dynamics 
generated by H, [33]. This in turn implies that the right-hand side of (l54l) should be x°-independent. 
In order to see this, we substitute (1351) and fl52D in dSH) and use (l30l) and fl3T]) to derive 



((A^O) = 7^ HA(xO),D-i/2[^2^(j^.S)2]A'(x°)^ 
2a^ 

^ (A(x°) + A°(x°)), Z}-i/2[^ _ . s)2](A'(x°) + A'°(x°)) >- 



K 

2m 



^A(x°),D-i/2e'(x°)^ - ^E(x°),D-i/2A'(x°)^ , (55) 



where E{x^) = — A(x°) — i^A^{x^). We can use ( ITTj) or equivalently (l30ll to check that the x 







derivative of the right-hand side of (l55l) vanishes identically. Therefore, (l55l) provides a well-defined 
inner product on V. Endowing V with this inner product and (Cauchy) completing the resulting 
inner product space we obtain a separable Hilbert space which we shall denote by Ti. This is the 
physical Hilbert space of the relativistic quantum mechanics of the Proca fields. 

The inner product (l55l) is identical with an inner product obtained in [29]. We will show in 
Section 4 that it is a special example of a larger class of invariant inner products and that it has the 
following appealing properties: 

1. It is not only positive-definite but relativistically invariant §| 

2. Its restriction to the subspace of positive-frequency Proca fields coincides with the restriction 
of the indefinite Proca inner product fj^Tl) to this subspace. 



As seen from f l5^ . the operator f/j,o for any value of x° G M is a unitary operator mapping 7i to 
/C. Following [26t 127] we can use this unitary operator to define a Hamiltonian operator h acting in 
Ti. that is unitary-equivalent to H. Let Xq G M be an arbitrary initial x^, and h : 7i ^ 7i be defined 
by 

h:=U;o^HU,o. (56) 
Then, using (l29ll . (135|1 - (138|1 and (l53l) . we can easily show that for any A G V, 

hA = ihA, (57) 



^We give a manifestly covariant expression for this inner product in Section 5. 
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where A is the element of V defined by: 74(x°) := A{x'^), for all x° G M. As discussed in [22] for 
the case of KG fields, one must not confuse (1571) with a time-dependent Schrodinger equation giving 
the x''- dependence of A{x^). This equation actually defines the action of the operator h on the field 
A. h generates a time-evolution, through the Schrodinger equation 

ih-^A^o = hA^o, (58) 

that coincides with a time-translation in the space V of the Proca fields; if Aq = A^o is the ini- 
tial value for the one-parameter family of elements A^o of V, then for all G M, A^o(x^') = 
(^e~^^^°~^o)^^^Ao){x'^') = Aq{x^' +x^ — Xq). Furthermore, using the fact that U^o is a unitary operator 
and that H is Hermitian with respect to the inner product fl52|) of /C, we can infer that h is Hermitian 
with respect to the inner product (IMl) of Ti. This shows that time-translations correspond to unitary 
transformations of the physical Hilbert space Ti. 
Next, we define U : H H' a.nd H' : H' ^ W by 

U:=pU,o (59) 
H' ■=UhU-^ = pHp-\ (60) 

where p is the unique positive square root of rj^, i.e., p := It is not difficult to see that 

P = 7-—= , (61) 




2m^ y P_ p+ J 2a4^ 

where p± : Ti ^ Ti are Hermitian operators given by 

Pi := {^^M T 1) [D^'^ - mD~^'^] i)^ + ^M^D-^I^ ± D^l\ (62) 

and \) is the helicity operator (12^ . We can check that the operator p viewed as mapping /C to H! is 
a unitary transformation; using p^ = p = ^/r]^, we have (p^,pC) = C ^ ^- This in 

turn implies that U : 7i ^ Ti' is also a unitary transformation, 

{UA, UA') = {{A, A')) , VA, A' en, (63) 

and that H' must be a Hermitian Hamiltonian operator acting in Ti,'. 

We can compute H' by substituting (EI]), (I62l) and (IMl) - (EH]) in (EOl). This yields 




H' = h [ ^ ] = hVDEs, (64) 



which is manifestly Hermitian with respect to the inner product (■, ■) of 7i'. The Hamiltonian H' 
is precisely the Foldy Hamiltonian pTl [16] . Here we obtained it by a systematic application of the 
methods of PHQM [5]. 



13 



Next, we derive the explicit form of U and its inverse. Using fl59|) . fl6Tl) . fl62l) . fl53l) . and fl35l) . we 
have for all A G 7^, R 



UA=- f^i ^A(xO)-zil-iE(xO) 
2 V V HA(xO) +zil-iE(x°) 



(65) 



where the operator ii: Ti ^ Ti. and its inverse are given bj0 

H = [DV4 _ ^/)-l/4][^2 ^ _^^-l/4^ ^-1 ^ ^JJ-I/A _ _^-l/^l/4][^2 ^ _^-l^l/4^ (qq) 

The inverse oiU is also easy to calculate. Let ^ G 7i' be a six-component vector (as in (15T|) ) with 
components and Then in view of 0651) . lA^^C, is the Proca field A G satisfying the following 
initial conditions: 



H4) 



D-1/4 



ni/4 



H4) = 



IM 



(67) 
(68) 



where we have made use of equations ([15]), ((301), ([33]), ([35]), ([53]), ([HI]), and ([62]). By virtue of p71), 
for all x° G M, we have 



A(x°) 



'1 „-i(xO-xO)DV2^i ^i(xO-zO)DV2^2' 



(69) 
(70) 



In complete analogy with the case of KG fields [27J, we find that the pairs (7Y, h), (/C, if), and (7-^', iJ') 
are unitarily equivalent. Therefore they represent the same quantum system. 



3 PT, C and CPT-symmetries of Proca Fields 

According to [35j, the generalized parity (P), time-reversal (T), and charge grading or chirality (C) 
operators for the six-component Proca fields are given by 



3 

^:=EE / ^'ka,,.(k)|<|.,,.(k))(<|.,,.(k)|, 

3 

^ := EE / a,,,(k) |$,,^(k)) ^ (<|.,,.(k) 

3 

C:=EE / f^'kor,,<,(k)|^,,^(k))($,,.(k)|. 



(71) 
(72) 
(73) 



-"^^Transformation is known as the Foldy's transformation, although it differs shghtly from the expression given 
by Foldy fT^ and Case [T^, namely Wpoidy ^3^- This difference turns out not to have any effect on the definition 
of the physical observables such as the position operator. 

"'^"'^It is interesting to see that the arbitrary parameter 7, introduced in the six-component formulation of the Proca 
equation, does not appear in (jM]) and ([65]). 
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where cre.o-(k) are arbitrary signs (±), and ^ is the complex-conjugation operator defined, for all 
complex numbers z and state vectors ^P, $, by i<z := z* and (^(<I>|)|\E') := ^($|\E') = (^|$). Note 
that we use the biorthonormal basis {\E'e,cr(k), '^'e,cr(k)} to define P, T, and C. This basis can be 
obtained by replacing the circular polarization vectors Me,/i(k) by the linear polarization vectors 
a,(k,cT) in (jli]) and (|i6|). 

Clearly there is an infinity of choices for the signs cTe^o-(k) each leading to a different P, T, and C. 
Following [27] we will adopt the natural choice associated with the label e appearing in the expression 
for the eigenvalues and eigenvectors of H . This yields 

^:=E E / rf'k6|<|.,,,(k))($,,,(k)|, (74) 
^:=E E / c?^ke|<|.,,,(k)) ^ ($,,_,(k)|, (75) 
C:=$^ 5^ / ci3ke|vl/,,,(k))(<|.,,,(k)|. (76) 

The consequences of these relations are identical with those obtained in [27J for the KG fields. Because 
of the close analogy with the case of KG fields, here we omit the details and give a summary of the 
relevant results. 

Substituting (il) and (iEl) in (CI) - ([76]) yields 

P = S3, T=S3*, Vr = ^. (77) 

Hence, the PT-symmetry of the Hamiltonian (l36ll means that it is a real operator. Note that 
= = So and VT = P ■ T, which is a direct consequence of ( 148]) , [35] . Similarly, we find 

C = ^( = h~'D-'l'H = (78) 

" ' -H2 -Hi / v^' 



which in view of (1491) is consistent with the identity C = ri^^V (equivalently 77+ = VC o) [35j . 
According to (178]) . C is a Z2-grading operator for the Hilbert space that splits it into the spans of 
the eigenvectors of H with positive and negative eigenvalues, respectively. 

We can use the unitary operator p : K. ^ H' to express the generalized parity, time-reversal, and 
chirality operators in the Foldy representation, namely V := pVp~^, T := pT p~^, and C := pCp~^. 
The symmetry generators VT' and C that commute with the Foldy Hamiltonian H' have the form 

V'T = VT = i., C = = S3. (79) 

Clearly, the P'T'-symmetry of H' is related to the fact that H' is a real operator, and the C'-symmetry 
of H' arises because it is proportional to C. 



^^Note that V, T and the metric operator 77+ depend on the choice of the biorthonormal system, while C is 
independent. In ^45J this fact has not been considered in the factorization of the metric operator. 
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Similarly we can use the unitary operator U : Ti ^ Ti' to define the generalized parity, time- 
reversal, and chirality operators for the ordinary Proca fields A ETi. This yields 

P := U-^V'U, T := U-^T'U, C := U'^CU. (80) 

Using Eqs. (1^3]) - fl68l) . and fll7p . we may obtain explicit expressions for P and T. Here we give the 
corresponding expressions for PT and C, that actually generate symmetries of the Hamiltonian h. 
The result is, for all A G 7i and G M, 

[PTA](a;°) = (-A°(-x°)*,A(-x°)*), (81) 
[CA] (x°) = «D-^/2^(x°) =: Ae(x°). (82) 

These relations imply that, similarly to the case of KG field [27], PT is just the ordinary time- 
reversal operator [1^, and "t/ie PT -symmetry of h means that the order in which one performs a 
time-translation and a time-reversal transformation on a Proca field is not important." Furthermore, 
chirality operator acting in Ti is simply given by C = iL>-i/2^o, which is a Lorentz scalar [6] . Recalling 
C^A = A, we observe that C : V ^ V is an involution, = 1. Therefore, we can use it to split V 
into the subspaces V± of ±-energy Proca fields according to V± := {A± G V|CA± = ±A-|-}. Clearly, 
for any A G V, we can use C to introduce the corresponding i-energy components: 

A±:=^(A±CA)GV±. (83) 

Clearly A = A+ + A_, so V = V+ © V_. Restricting the inner product fl55l) to V± (and Cauchy 
completing the resulting inner product spaces) we obtain Hilbert subspaces H± of H. 

Next we recall that C is a Hermitian operator acting in Ti.'. Thus, in view of ( !80l) . and ( l63l) . we 
have {{A, CA')) = ([CA, A')). Therefore, C : 7^ ^ 7^ is a Hermitian involutioill!] and for all A± e V±, 
((y4+. A-)) = 0. This in turns implies 7i = 7Y+ ©7Y_. The generalized chirality operator C is actually 
the grading operator associated with this orthogonal direct sum decomposition of Ti. In other words, 
similarly to the case of KG fields [27], C is a Hermitian involution that decomposes the Hilbert space 
into its ±-energy subspaces. As a result, "t/ie C-symmetry ofh means that the energy of a free Proca 
field does not change sign under time-translations." 

4 The Most General Admissible Inner Product 

The metric operator rj^ and the corresponding invariant positive-definite inner product depend on the 
choice of the biorthonormal system {^^ /^(k), $e /^(k)}. The most general invariant positive-definite 
inner product corresponds, therefore, to the most general biorthonormal system that consists of 
the eigenvectors of H and W . In the construction of {\E'e,/i(k), $e,/i(k)} we have chosen a set of 
eigenvectors of H that also diagonalize the helicity operator A. As a result, both H and A are 
77+-pseudo-Hermitian. This in turn is equivalent to the condition that H and A are among the 

^•^This marks its similarity to the chirality operator for spin 1/2 fields. 
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physical observables of the system. In the following we construct the most general Lorentz invariant 
positive-definite inner product that fulfills this condition!^ 

The most general positive-definite operator fj_^. that renders both H and A r/_|_-pseudo-Hermitian 
has the form fj^ := A^rj^A where A is an invertible linear operator commuting with H and A [39|[35]. 
We can express A as 

^ = Y.ll I ^^'ka,,,(k)|*,,,(k))(<l>,.,,(k)|, (84) 
where a^^hi}^ are arbitrary nonzero complex numbers!^ This implies 

7^+:=^W = 5^ I ^^'k|a,,,(k)n$,,,(k))($,,,(k)|. (85) 

Substituting (jH]) and fH6l) in fl85|) . and carrying out the necessary calculations, we find 

?7+ = [L+M + {Lt - L^_)as - L^N] f)^ + [L;Af + LZcis] f) + L^N + ag, (86) 



where 



LI := i| ci3k[|a+,i(k)p + e|a+,_i(k)p±|a_,,(k)p±e|a_,_i(k)p] |k)(k|, 
Ll := 1| d3k[|a+,o(k)P±|a_,o(k)P] |k)(k|, 



^7) 



( 72^ + 1 ^^D - 1 
27 1^ - 1 + 1 

D-V^ / 72^4^^ 72^4 _^ 

27A^2 1 72_^4_^ 72^4^ + /} 



I in . , . I' (89) 



Having computed 77+ we can easily obtain the corresponding inner product according to 
77+-). This yields, for all ^, C G 7i', 



UX))n^ = \ [7^^+,e+,+iC+^ +7"'-<^-,e+,_iC_H + ^^'>0-,oCV - ^^^e_,oCV, (91) 

where the operators 6e,/i ^H. are defined by 

e,,^ := [L+D'^/^ _ ^2/^^o^-V2][^2 ^ 2.7Z}^/2^ + A^2/^^o^-V2_ (92) 

If we view Ti' as a complex vector space and endow it with this inner product, we obtain a new inner 
product space whose Cauchy completion yields a Hilbert space which we denote by /C. Finally, we 
recall that, defining ((A,A')){(i} := ((?7a.oy4, f/j,oA'))^^, we can endow the space V of the solutions of 
Proca equation with the positive-definite inner product 

((A,A'))m := ^ [^A(x°),0+,+iA'(x°)^ + ^E(x°),0+,_iE'(x°)^ 

+z{^E(x°),e_,oA'(x°)^ - ^A(x°),0_,oE'(x°)^}] , (93) 



^^This amounts to restricting the choice of the metric operator in the spirit of [51]. It does not however fix the 
metric operator (up to a trivial scaHng) because H and A do not form an irreducible set of operators. 

i^The condition that H is 77+-pseudo-Hermitian yields A = Y.e=± Y.h=o,±i I]/i'=o,±i /rs '^^^'^h,h' {^)\^<i,h (k)) {<^e,h' (k) | 
which does not commute with A unless ^, — Sh.h'Cte,h- 
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which is dynamically invariant. This fact may be checked directly by differentiating the right-hand 
side of ( |93l) with respect to x° and making use of (fTTil or equivalently (!30l) . 

By construction, fl93l) gives the most general positive-definite inner product that renders the gen- 
erator of time-translations as well as the helicity operator Hermitian. We next impose the condition 
that (I93p be Lorentz-invariant as well. This will give rise to the most general physically admissible 
inner product on the space of Proca fields. 

To quantify the requirement that the right-hand side of (!93|) be a Lorentz scalar we compute the 
inner product of two plane- wave solutions of the Proca equation, i.e., (127|) . Using fl2^ . (!32l) . (!87|) . 
(!88ll . and ( j92l) in ( j93l) . we find after a rather lengthy calculation 

((A,,,(k),A,,,,,(k'))){a} = -\N,,hm'a,,hu;,,S%k-k')6,,,,6h,h', (94) 

M 

where we have introduced the abbreviation a^ ^ := la^^hP- Because a^^h are positive real numbers and 
the right-hand side of flMl) does not involve x", this equation provides an explicit demonstration of 
the invariance and positive-definiteness of the inner product ( l93l) . 

Next, we recall from Refs. [3l HE] that relativistically invariant normalization of two plane- wave 
solutions is given by: ([A^^h(k; x), A^i^h'(k'; x))) = 2c<Jk5^(k — k') 6i,^^'6h,h'- This together with the 
fact that N^^h{k) must be a Lorentz scalar so that A^^h{k;x) is a four-vector, shows that ae,h are 
dimensionless positive real numbers, i.e., they do not depend on k and obey the same Lorentz 
transformation rule as scalars. This result ensures the relativistic invariance of the inner product of 
any two solutions A and A' having the general form: 

A = Y,J2[ d'kc,,h{k)A,h{k), (95) 

e=± h=0,±l -^^^ 

where c^^h{k) are complex coefficients. In view of flMl) and fl95l) and the fact that the inner product 
fl93|) is a Hermitian sesquilinear form, we obtain 

((^'^'))W = ^E E / rf'ka,,^k|iV,,(k)|^<,(k)<,(k). (96) 

Because A and A^^h{k) are four- vectors and d^k/u^ is a relativistically invariant measure [3l US], 
Ce^h{k) obeys the same Lorentz transformation rule as co^'^- This in turn implies that in order for the 
inner product (jHSD to be scalars, a^^h must transform as scalars. Therefore, in view of (|57|1 and (155]) . 
we conclude that and are dimensionless real numbers given by 

Inserting these relations in ( !92l) and using ( !30l) we find 

0+,+iA = 0+,o/^-i/2E, 0+,_iE = -e+,oi^-'/'A. 

These in turn allow us to express the most general physically admissible inner product on the space 
V of the Proca fields, namely (jHSD, as 

((A,/1'))m := ^ [-<A{x'),e+,oD-'/'E\x')y- ^E(x°), 0+,o Z^-^/^A'(a;°) ^ 

-I A(x°), e_,oE'(a;°) ^ - ^E(x°), e_,oA'(x°) ^}] , (98) 
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where, according to Eq. fl92l) . 

e±,o = [Lt - L° + L^f) + L° . (99) 

Note that is a positive number and we can absorb it in the definition of the field. Therefore, 
the inner product flUHl) actually involves five nontrivial free parameters. Endowing V with this inner 
product and Cauchy completing the resulting inner product space we obtain a separable Hilbert 
space which we denote by 7i{a}- 

We can express the inner product (pH!) in terms of the indefinite Proca inner product (HTI) . First, 
we use the unitary operator U^o : H{a} ^ /C to find an expression for the action of the helicity 
operator 9j := U'^hU^o on a Proca field A G 'H{a}- [^^{x^) is a Proca field whose components can 
be evaluated using (l53l) . fl33|) . and f|T7|) . This gives 

[^lA](x°)= (0,[)A(a;°)). (100) 

Because [A, iJ] = and and h are respectively related via a similarity transformation to A and 
H, we have [i3,h] = 0. This together with the identity S^^ = suggest to use 9) to split V into the 
subspaces V±i and Vq of ±1- and 0-helicity Proca fields according to Vh '■= {Ah G V|i3A/i = hAh} 
where h G { — 1,0, +1}. We can use to define the ±1- and 0-helicity components of the Proca fields 
A eV according to: 

A±i := i {Sj^A ± 9jA) G V±i, Ao := {A - 9)^ A) G Vq. (101) 

Clearly, A = A^i + A^i + Aq which implies V = V+i © V_i © Vq. Restricting the inner product 
( 198|) to V±i and Vq (and Cauchy completing the resulting inner product spaces) we obtain Hilbert 
subspaces 'H{a}±i and Hiajo of TYjn}. 

Next we recall that is obtained via a unitary similarity transformation from the Hermitian 
operator A. This in turn implies 

iiA,f)A')){a} = im,A')){a}. (102) 

Therefore, ^ is a Hermitian operator acting is 7i{a}- It is a physical observable that measures the 
helicity of A. Furthermore, in light of (11021) and (1 10 II) . 

((A+i, A_i)){„| = ((A+i, Ao)){a} = {{A-l, Ao)){a} = 0, VA eVh, he {-1, 0, +1}. 

These relations show that Hia} = "^{aj+i © '^{a}-i © "^{ajo- Also, as we have shown in Section [3l 
the chirality operator C is a Hermitian involution acting in ?i{a} which decomposes the Hilbert 
space 'H{a} into its ±-energy Hilbert subspaces Ti.{a}±- Thus, using the operators Sj and C, we can 
decompose the Hilbert space into six mutually orthogonal subspaces H{a}{t,h), where e = {+, — }, and 
hG {— 1,0,+1}. Employing this decomposition of A in ( l98l) . we find 

((A,A')){a} = 5^ ea,,UiA,H,KJ,^. (103) 

€=± h=0,±l 
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Therefore, fl98l) is a linear combination of the indefinite inner products fHTl) of the components A^^h of 
A with suitable non-negative coefficients. Note that the decomposition of the fields into the helicity 
components fllOll) is clearly frame-dependent and so (11031) hides the Lorentz covariance of the theory. 



But the inner product (pHjl does not involve the explicit splitting of the fields into (e, /i)-components. 

In the remainder of this section, we obtain the Foldy transformation associated with the general 
inner product fl^Sl) . First, consider the invertible linear operator A of This is a unitary operator 
mapping K. onto /C, for it satisfies, for all ^, C G /C, 

{{A^, AOK = (-^e, V+M) = (e, A^V+M) = Ch^- (104) 
We can find the explicit form of A and A^^ by substituting flH|) and fHB]) in flM|) . This gives 



A = [F^asM - FVgAT + (F+ - F°)(To][)' + [Fla^M + F-ao]i) + F^a^N + F°ao, 
A-' = [F^a^M - F\zN + (F+ - F°)ao]f)' + [Fr^gM + p-a^]^ + FVaiV + F°ao, 



where 



2 L 

1. 



Next, we define p : JC ^ H' and U{a} '■ '^{a} — ^ by 



In view of = p = and fl85l) . p is a unitary transformation mapping /C onto 7i', 

(p^,pC) = (e,pV~c) = (e,^Wc) = ((e,c))^+- 

This in turn implies that U^a} '■ ^{a} ^ is also a unitary transformation, 

{U^,}A,U{,}A') = {{A, A')){„}, VA, A' e 



(105) 
(106) 



(107) 
(108) 
(109) 



(110) 
fill) 



(112) 



(113) 



We can compute the explicit form of the unitary operator U{a} and its inverse. This requires 
computing p and its inverse. Substituting ( l6Ti) . ( l62l) . f llOSp and f ll06p in (11 101) . we have 



P 



P+,+ P-,+ 



V p™,- p^ 



p+: 



-p- 



-pj°; P+°! 



(114) 



where the entries are the following operators that act in Ti. 



Pe 



Pe,e' •" 



(115) 
(116) 
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Next, we use ffTTTD . ffTTD . ffTT5|) . (1531). and ([35]), to find, for all AeTi, 



^^"^"^ 2 V ^ \^ ii_,+ A(xO) + zil_,_E(xO) 
where the operators 11^ e' : Til ^ Til, are given by 



(117) 



'118) 



Again the arbitrary parameter 7 drops out of the calculations. = U^J^p^^ is also easy to 

calculate. Let ^ G 7Y' be a six-component vector (as in (15T|) ) with components and Then in 
view of (11171) . is the Proca field A G Ti-ia} satisfying the following initial conditions: 



r)-l/4 



'MK 



"+,0 "-,0 



A°(a;°) = -z 



"-,0 



(119) 



A(x°) = ./^ (y;^^^^ + HI^^C^) , A{4) = -r Md^^ {^-\e - ^--\e) , (120) 
where is the inverse of 11^^+ given by 



iiri 



/)l/4 



f)2 + ZlD-l/^f^ + 



DV4 



:i2i) 



and we have made use of (USD, ([3SD, dSSj), dSSj), and (fTTlD - (fTTBD . By virtue of ([IT]), for all 



Xq G M, we have 



A°(x°) 
A(x°) 



D-V4 

\/MK 
M 

K L 



(122) 
(123) 



Next we recall that since [-ff, .4.] = 0, we have h{n} = h. Therefore, in hght of the above analysis, 
the pairs (7i{a},h), {IC,H), and (Ti.',H') are unitarily equivalent; they represent the same quantum 
system. 

Finally, we give an alternative form of the Foldy transformation (I117p . In view of (|83|) and 



K I ill ,A,fxS) 



il-lE = we have 



In particular, setting a^./i = 1, for all e G {+, — } and h G { — 1, 0, +1}, yields 



(124) 



UA 



il 



A+(x[]) 



M \ A_(Xo) 



(125) 



where A± is i-energy component of Proca field. These equations show that ile^+A^ (specially ilA^ 
satisfy the Foldy equation [T71 [T6] 



(126) 



21 



5 Conserved Current Density 



The relativistic and time-translation invariance of the positive-definite inner products (|98|) suggest 
the existence of an associated conserved four- vector current density J'^^y In this section we use the 
approach of [6j to compute J'^^-^ for the case ^ = 1 and obtain a manifestly covariant expression for 
•^fa} with ae,h = 1 which for brevity we denote by J'^. 
According to (jHSD, for all A e V, 

((A^)){a} = ^^d3x{A(xO,x)*.(e+,oI)"'/'E(xO,x)) -E(xO,x)*.(e+,oI?-V2A(xO,x)) 

-i[A(xO,x)*.(G_,oE(a;0,x)) -E(x°,x)*.(e„,oA(x°,x))]}. (127) 

In analogy with non-relativistic QM, we define the 0-component of the current density J^^y associated 
with A as the integrand in (11271) . i.e., 

4a}(^) ■■= ^ {MxT ■ {e^,oD-'/'M^)) - E(a;)* ■ (e+,o/^-i/^A(x)) 

-z [AixY ■ (e_,oE(x)) - E(x)* ■ (0_,oA(x))]} . (128) 

Setting a^^h = 1 for all e G {+, — } and h G { — 1, 0, +1}, yields 

jO(a;) := Jl_ \A{xy ■ D-'^^Eix) - E(x)* • D-'^^A{x)} . (129) 
2m l J 

In order to obtain the spatial components of J'^, we follow the procedure outlined in Ref. [6]. Namely, 
we perform an infinitesimal Lorentz boost transformation that changes the reference frame to the 
one moving with a velocity v: 



/O 

X —>■ X 



a;° - /3 ■ X, X ^ x' = X - (3x°, (130) 



where f3 := v/c. Assuming that J'^ is indeed a four- vector field and neglecting the second and higher 
order terms in powers of the components of f3, we then find 

J\x) J'\x') = J\x) - (3 ■ 3{x). (131) 

Next, we use (11291) to obtain 

:= — |a'(x')* ■ D'-'/^E'ix') - E'(x')* ■ D'-'^^A'ix')] , (132) 

where x' := (x'°,x'), D' = -V'^+m'^, E'(x') = -A'{x')-VA'\x'), and A'(x') := dA'{x')/dx'^ . This 
reduces the determination of J to expressing the right-hand side of fll32p in terms of the quantities 
associated with the original (unprimed) frame and comparing the resulting expression with (11311) . 
Under the transformation (11301) . the four- vector A = [A^,A) transforms as 

A°(x)^A'°(x') = A°(x)-/3-A(x), A{x) ^ A'{x') = A{x) - (BA'^ix). (133) 

Therefore, in view of (I130p . we can easily obtain the transformation rule for E. The result is 

E(x) E'(a;') = E{x) + /3 x (V x A(x)). (134) 
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Because the chirality operator C = iD ^^"^do is Lorentz invariant [6], 

D'-'/^d'o = D-'/^do. (135) 
Substituting fll33p - f ll35p in f ll32p . making use of f ll3ip . we obtain 

r{x) = ^ \a*{x) ■ d'D-^'^A{x) - [d'A{xy] ■ d-^'^a{x)- 

2m I- 

[A*{x) ■ d]D-^/^A'{x) + [D-^/^A{x) ■ d]A\xy} , (136) 

where d := {d^, V) and for any two four-vectors Vi and V2, f 1 ■ V2 '■= v^V2^. This relation suggests 

J^^{^x) = — {a*{x) ■ d^D-'^^A{x) - [d^'Aix]*] ■ D-^/^A{x)- 

[A*{x) ■ 9]D-i/'>(x) + [D-'^^A{x) ■ 9]A^(a;)*} . (137) 

It is not difficult to check (using the Proca equation) that the expression for J° obtained using this 
equation agrees with the one given in 



We can use ([2]) and (1521) to further simplify fll37p . This yields 

■J^i^) = ^ {M^TF^^i.^) - F-^{xrA,,{x)} , (138) 
zm 

where F^^ := d^A'^^ — A^^. The current density which is generally complex- valued has the 
following remarkable properties: 

1. The expression (11380 for is manifestly covariant; since A and A^ are four- vector fields, so is 

2. Using the fact that both A and A^ satisfy the Proca equation (jl]), one easily checks that 
satisfies the following continuity equation. 

d^J^" = 0. (139) 

Hence it is a conserved current density. 

Next, we use fll38p to derive a manifestly covariant expression for the inner product flSSp on the 
space of solutions of the Proca equation (ITTi) . The result is 

{{A, A')) = daix) n,ix) {^.(a:)*^^ (a:) - F'^^(x)M^,(x)} , (140) 

where a is an arbitrary spacelike (Cauchy) hypersurface of the Minkowski space with volume element 
da and unit (future) timelike normal four- vector n^. Note that in deriving fll40p we have also made 
an implicit use of the fact that any inner product is uniquely determined by the corresponding norm 

nzi. 

Using the same approach we have calculated J^^^j for a^^h 7^ 1 and checked that indeed it is a 
conserved complex-valued four-vector field. But we were not able to obtain a manifestly covariant 
expression for J'^^y in this case. As the expression for J'^^y is highly complicated we do not present it 
here. 
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6 Physical Observables and Wave Functions for Proca Fields 



The unitary equivalence of the representations (7i{o},h), {IC,H), and {H',H') for the quantum sys- 
tem describing the Proca fields allows for the construction of the observables of this system using 
any of these representations. Because H' = n®n andH = L'^(R^) © L'^(R^) © L'^iR^), the Foldy 
representation {7i', H') is more convenient for this purpose. In this section we construct the observ- 
ables in this representation and use the unitary map tl{a} '■ ^{o} ~^ to obtain their form in the 
standard (covariant) representation (7Y{a},h). Again we follow closely the approach used in [27j to 
construct the observables for KG fields. 

In the Foldy representation, we introduce the following set of basic observables 

X'^:=x©S„, P'„:=p©S„„ 5^ := 1 © S„„ (141) 

where, x, p = /l.^, and 1 are respectively the position, momentum, and identity operators acting in 
L^(]R^), m G {0, 1, 2, ■ ■ ■ , 35}, and S^'s are given in Eq. ([9]). In the following, we will omit '1©' for 
brevity. In particular, we will identify S'^ with S^. 

As we mentioned in Section 1, the operators Xq, S3, S12, S15, S32, and S35 form a maximal 
commuting set of observables acting in Ti' . In particular, we can use their common eigenvectors, 
namely 

e(^'^):=|x)©e,,„ xGM^ eG {-,+}, sG{-1,0,+1}, (142) 

to construct a basis of Ti' . In ( 11421) . e^^s are the vectors defined in (ITUl) and |x) are the 5-function 
normalized position kets satisfying x|x) = x|x), (x|x') = 5'^(x — x'), and Jj^g (i^x |x)(x| = 1. It is 
easy to see that indeed Xq^x'*^= x^x''^-*, S3 e^x and Si2^x s^^'^\ Furthermore, 

(ei^'^US'^'^) = ^..'5s,.'<5^(x-x'), EE / ^'xi4^'^^)(ei^'^^i = So. (143) 

We can express any six-component vector ^' G Ti' in the basis {^x'*''} according to 

^' = E E / rf'x/(e,.,x)e4^'^), (144) 

where / : {— ,+} x { — 1,0, +1} x R-^ ^ C is the wave function associated with ^' in the position- 
representation, i.e., 

/(e,s,x):=(ei^'^),v[/'). (145) 
Also, the action of a physical observable O' : Ti' —>■ Ti' on the state vector G Ti' is given by 

^'^' = E E / rf'x[6/(e,s,x)]^(^'^), (146) 

where 

6/(e,s,x):=E E / ^'x'/(e',.',x')(4^'^),0'ejyV (147) 
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This provides the representation of observables in terms of hnear operators acting on the square- 
integrable wave functions /. 

Next, we introduce the operators 

Xm{a} := ^"jX'^^a}, Pm{a} := P'm^a}, Sm{a} ■= ^{a}^mi^{a} , (148) 

that act in Ti.{a}, and define the Proca fields 

4aU = (^?&- ^Wx) - ^{-.1 (149) 

which form a complete orthonormal basis of Ti^a}'- 

((4a}t4a}x'))w=^-'^-'^'(---')' EE / ^'-l4a}t)(4a}tl = -0{a}. (150) 

Here we have used (I113p . fll43p . and (11491) . So{a} coincides with the identity operator for H{a}, and 
for all A, A' e n{a}, the operator \A){A'\ is defined by \A){A'\A" := {{A', A")){a}A, for any A" e 7^{a}. 
In view of (11501) and 

Utll ^)){"} = i^wAall^i^}^) = (^x '^^ = /(^' ^' (151) 

we can express any Proca field A G 7i{a} in the basis according to 

^ = E E / rf'x/(e,.,x)4„t (152) 

Here we do not label the wave functions / with the subscript {a}, because they do not depend on 
the choice of the parameter a^^h- The proof of this assertion uses the unitary operator 

W^-^JW, (153) 

that maps Ti. onto "^{0} and is identical with the one given in [6] for the KG fields. 

The physical observables 0{a} '■ '^{a} ^{a} are uniquely specified in terms of their representation 
in the basis {^j^jx}' ^ ^ '^{c} 

^{a}^ = EE / ^'x[6/(e,.,x)]4„j)^, (154) 

where 0/(e,s,x) is defined by ([T47D. This follows from {{A^(;^l,O{a}A['^0){a} = (d' 0'^^,' ''^ 
In view of (11501) and (11521) . we can express the transition amplitudes between two states (the inner 
product of two Proca fields) in the form 

{{A,A')){,y = Y, I d'^f{e,s,^rfie,s,^). (155) 

More generally, for any observable 0{a} '■ T^la} "^{a}) we have 

((Ao{a}A')){a} = E E / c?^x/(e,s,x)*6/'(e,s,x) 

e=± s=0,±l "^K' 

= EE / rf'x[0/(e,.,x)]7'(e,s,x). (156) 



e=± s=0,±l ■ 
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The above discussion shows that we may view the wave functions / as elements of and similarly 
to the case of KG fields , formulate the QM of Proca fields in terms of these wave functions. 

In this formulation the observables are Hermitian operators O acting on the wave functions. For 
example, the action of xo{a}, Po{a}) S3{a}y and Si2{a} on A corresponds to the action of the operators 
xo{o}, Po{a}, ss{a}, and Si2{a} ou /, where 

Xo{a}/(e, ■5,x) := x/(e, s,x), po{a}/(e, s, x) := -inV/(e, s, x), (157) 

■S3{a}/(e, s,x) := e/(e, s,x), Si2{a}/(e, s, x) := s/(e, s, x). (158) 

These equations show that /(e, s,x) are the position wave functions with definite chirality (sign of 
the energy) e and spin s (say along the x^-direction). They are however not the eigenfunctions of the 
helicity operator. To see this, we recall from fHOl) that 7i' = 7i © 7i, and the spin operator acting in 
Ti. is given by S of Eq. fl25|) . This shows that the spin operator acting in Ti.' io 

S':=(E28,-S2o,S8)= s)' ^^^^^ 

Denoting the spin operator acting in H^a} by Sj^}, we have 

S{„}=W^-;jS'Wm. (160) 
The projection of this operator along the momentum po{a} gives the helicity operator: 

(P0{a} ■ S{a})/(e, S\ X) = Eijk Oj /(e, s^, x), (161) 

where {s^, s^, s^) = (+1, —1,0), and we employed fll47p . 

Similarly, the action of the Hamiltonian h on A corresponds to the action of the operator h : = 



■53{a} yPojo} + "^^c^ on the wave function /: 

h/(e, s, x) = heV-V^ + M^f{e, s, x). (162) 
Subsequently, in view of (]162p . and (|58l) . the dynamics of the evolving Proca field A^o is determined 

{a}: 



in terms of the wave functions /(e, s, x; x°) = {{A^s'^f^^, Axo))^a}, according to 



ihdofie, s, x; a;°) = eV-h^V^ + m^c^ /(e, s, x; x°). (163) 

Furthermore, applying ido to both sides of fll63p . we can check that the wave functions / also satisfies 
the KG equation: [d^ - + M^]f{e, s, x; x°) = 0. 

Next, recall that because the time-reversal operator ( !79l) acting in Ti' commutes with Xq, the 
eigenvectors ' may be taken to be real. In this case the action of the time-reversal operator 
T = PT on any A G ^H^a} is equivalent to the complex-conjugation of the associated wave-function 
/, i.e., T/(e, s, x) = /(e, s, x)*. Similarly, we can identify the operators S3 and S12, respectively, with 
the chirality and spin operators acting on the wave functions /. 

As we shall see in the following section the wave functions /(e, s,x) furnish a position represen- 
tation for the QM of the Proca fields. The corresponding position operator is the spin-1 analog of 
the Newton- Wigner position operator for the KG field |211 [22] and similarly to the latter fails to be 
relativistically covariant,^ This means that the above-mentioned position-representation provides 



^^Throughout this paper we express spin and angular momentum operators in units of h. 
^''See however HSl. 
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a non-covariant description of a quantum system that also admits a unitary-equivalent covariant 
description in terms of the Hilbert space 7i{a} and the Hamiltonian h. 

7 Position, Spin, and Localized States 

As discussed in [27] for the KG fields, the canonical quantization scheme that provides the physical 
meaning of the observables yields the Foldy representation of the quantum systemF^ This suggests 
that the operators xo{a} and Po{o} that clearly satisfy the canonical commutation relations may be 
identified with the position and momentum operators acting on the wave functions /. This in turn 
means that the operators Xq and Pq in the Foldy representation and the operators xo{a} and Po{a} 
in the (Tijn}, h)-representation also describe the position and momentum observables. In particular, 
the basis vectors ' and ^{n}x determine the states of the system with a definite position value x; 
they are localized in space. They also have definite charge or chirality (sign of the energy) and spin 
(say along the x^-direction) . 

By construction, v4|^J|^^ are delta-function normalized position eigenvectors, i.e., they satisfy (11501) 
and 

m^AZ = ^4^ (164) 

Similarly, we identify the chirality operator and the spin operator along the x^-direction with C = 
W|~|S3W{(j} and su = W|~|Si2W{ii}, respectively. This implies 

In view of Eqs. (I150p . (I164p . and (I165p . the state vectors ^j^i^ represent spatially localized Proca 
fields with definite chirality e and spin s. We can associate each Proca field A G "^^{0} with a unique 
position wave function, namely /(e, s,x). As we explained in Section [6l we can use these wave 
functions to represent all the physical quantities associated with the Proca fields. We also emphasize 
that according to (11611) . /(e, s,x) are not the helicity eigenfunctions. Therefore, the state vectors 
^{o}x '^'^ ^'^^ have definite helicity. This is in complete accordance with the Heisenberg uncertainty 
principle: since A^^^^^^ are localized states ( with definite position ) they do not have definite momentum 
and hence definite helicity. 

Next, we recall that two Hilbert spaces 7i{a} for different choices of the parameters {a} are unitary- 
equivalent. We can use (11481) . (I149p . and the unitary operator (11170 to obtain the exphcit form of 
the localized Proca fields and the physical observables acting in 7i{a}. The resulting expressions for 
the position and spin operators are highly complicated. Therefore, in the following we shall only 
derive the explicit form of these operators in the covariant representation (7i,h). 

-"^^In the classical limit, each component of a Proca field corresponds to a classical free particle of energy E = 
±-\/|7^ + m^c^. Upon quantization E — s- ±fi,\/— + which signifies the relevance of the Foldy representation. 
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7.1 Explicit Form of the Localized States 



As suggested by (I155p . the wave functions /(e, s,x) belong to L^(M^). Moreover due to the imphcit 
dependence of A^'^^ on the appearing in the expression for U, /(e, s, x) depend on Xq. As in the 
case of KG fields this dependence becomes explicit, if we express /(e, s, x) in terms of A directly. To 
see this, we first substitute (EZD, dHBD and (fTI2]) in (fTl^ to obtain 

Al(^'^')(x')=e^ e-(-"--o)^^^'i^Mx), A^r'\x'>) = . il-^e-(-°--S)^^/^ |x) ® e,., (166) 

y/MK V ^ 

where (s^, s^, s^) = (+1, —1, 0). We then use this equation and fl55|) to compute the right-hand side 
of dlSTl). This yields 

fie,s\^) = /Z e-(^°-S)^^^^ [ilA,(a:°,x)]\ (167) 

V M 

where A'^ is the definite-chirality (definite-energy) component of A^^ with chirality e. Because A'^ and 
consequently ilAe satisfy the Foldy equation (11261) . we have e*'^^^°~^o^^^^^ilAe(x°, x) = ilAe(xg,x). 
Inserting this in (11671) . we find the following manifestly Xg-dependent expression for /(e, s*,x). 



/(e,s\x) 



[ilA,(a;°,: 



(168) 



Next, we use (11661) to compute the value of the localized Proca fields Ay'^ ■* at a spacetime point 
(xO,x): 

4^'^')(a;°,x) := (x|4^'^')(x°)). (169) 
Doing the necessary calculations, we obtain 



Ay^^'^^\x) = iesinO cos4> Ii, Ay^"' (x) = iesinO sincj) Ii, A^^"'^' (x) = iecos 9 Ii, (170) 



,0(e-l) 



0(e,0), 



\ V3 



A^'-'\x) 



( V2 



/ V3 



(171) 



where x := (x°,x), 6 and are the polar and azimuthal angles representing the direction of x — y, 

Vl := I2 + sin^ 6'cos^ I^, V2 '■= \ sin^ ^sin(20) /a, ^3 := | sin(26') cos0 /a, 
f 4 := I2 + sin^ 6 sin^ J3, t>5 := | sin(2^) sin 0/3, vq := I2 + cos^ 6' /3, 



and 



h 
h 
h 

Qi{e,k) :-- 



|x-y| 
/ . ^ 1 /""^ 



K 27r2 

M 1 



rfA; r]i(e, fc) <j ^'''li^''',,, ^'^ + 1^2(fc) {M~^[k^ + A^^/' - 1} L 



ft 27r2 



(i/c i7i(e, /c) 



|x-y| 

k sin(A;|x — y|) 

|x-y| 



?>n2{k)) {M-^[e + M^Y'^ - I] 



exp [— ze(x° — a;Q)\/P~+A?] 

[P + m2]1/4 



sin(fc|x — y|) cos(/i;|x — y|) 
/c|x — yP |x — yP 



(172) 
(173) 
(174) 
(175) 
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Figure 1: Plots of Ji, I2 and /s in terms of the radial distance |x — y|. The radial distance is scaled 
with the Compton wave length m~^. 

For x° = Xq, the integrals on the right-hand side (11721) - (11741) can be expressed in terms of the Bessel 
K-function (-ft'n), and the Hypergeometric function {pFq). The result, for both e = — 1 and 1, is 
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where z := |x — y| and F is the Gamma function. Fig. [T] gives the graphs of /i,/2, and I^. They 
involve a 5-function-like singularity at |x — y| = 0. This is a manifestation of the fact that Ay'^^ are 
localized at y. 

As seen from Eqs. (11701) and (I17ip . the localized Proca fields, unlike the localized KG fields [S], 
depend on the angles 9 and (direction of x — y). But as for the KG fields their position wave 
functions involve the delta functions; the position wave function f(e,s,yc){^', s', x') for Ax'^\xq) has the 
form (5e,e'(5s^s"^^(x — x'). 

Finally, we wish to emphasize that to the best of our knowledge the explicit form of the localized 
Proca fields A^'^^ have not been previously given. It is remarkable that we have obtained these 
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localized states without pursuing the axiomatic approach of Ref. [21j. The latter gives rise to the 
Bargmann- Wigner localized fields [2T|, [20] • 



7.2 Position, Spin, and Helicity Operators 

Using the unitary map U and its inverse, we can obtain the explicit form of the position operator 
xq that is defined to act on the Proca fields A ETi. Note that x •= is a three-component field 
whose components satisfy both ( fTTl) and ( fT3l) . Therefore it is uniquely determined in terms of the 
initial data (x(a:o), x{xq)) for some Xq G M. One can compute the latter using (11481) . (|^ . ([SZD, dSH]), 
and the identities: D = ^ + m^, 

f)Sf) = ^f), (179) 

[X, 0] = ^ S - ^ 0, [x, f)^] = ^ (OS + SO) - ^ [)^ (180) 

and [F(.^),x] = —i^^F{^, where F is a differentiable function. This yields 

/(x°) = 2)AO(xO) + k{xl), x'ixl) = (?) - f ) i°(^S) - ^ A(xO), (181) 

X(x°) = XA{xl), x{xl) = - ^) A(xS), (182) 

where 

X := X - i + - M-^D-^/'V + - MD-^r']t,S + ^Xm-^D^^"- - 1)SD. 

2u |j^| |J^| 

_ iJ^ _ zi^(J^-S)2 z[S(j^-S) + (J^-S)S] _ (Di/2-^)(j^xS) 

Next, we employ f|T7|) to express x(3^°) terms of fll8ip and (11821) . Simplify the resulting 
expression, we find, for all G M, 

X°(a;°) = (2) - z(a;° - xl)^D-'d^} A°(a;°) + -l^A(a;°), (185) 
X(x°) = {X - i(x° - a;;])J^D-^9o} A(s°). (186) 

In addition to being a Hermitian operator acting in the physical Hilbert space ?i, the position 
operator xq has the following notable properties. 



1. It respects the charge superselection rule [39], for it commutes with the chirality operator 
C = S3. This is easily seen by noting that xq and S3 are respectively obtained via a similarity 
transformation (11481) from Xq and 1S3 = S3, and that according to (I14ip . [Xq, S3] = 0. 

2. It has commuting components, and it commutes with the spin operator. This is because it is 
obtained via a similarity transformation (I148P from Xq which has commuting components and 
commutes with S'. 
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3. Its X -derivative gives the relativistic velocity operator; i.e., we have: 

V := ^ = 2 [h, xo] = lU-^ \H\ X'] U = , P° C. 

As we noticed in [28], this result means that the physical momentum is po C. 

4. It has the correct nonrelativistic limit: As c ^ oo, xq x. 

Similarly, we can obtain the explicit form of the spin operator Sq acting on A G 7Y. Again, noting 
that 5^ := SqA is a three-component field whose components satisfy fill I) and (1131) . we can determine 
y in terms of the initial data {y{xl),y{xl)). Using ([65]), ([67]), ([68]), ffMHD . ffTTOD . and 

[S, f)] = ^ (-^ X S), [S, f)2] = ^ [(-^ X 8)0 + f)(i^ X S)] , (187) 

these initial values are given by 

^0(a;0) = M-'D-'/^R X A(xO), ^^(xO) = -.M-^Z^i/^il x A(xO), 

^(xO) = SoA(xO), J^(xO) = SoA(x[J), ^^^^^ 

where 

§0 = S + ^[1 - mD-'/'MK X S) + ^[-M-^/^^/^ - 1](^ X S)i) 

D + M^^ (DV2 _ ^) . S) z{(j^-S)(J^x S) + (j^x S)(j^-S)} 



(189) 



2^^1/2 2A4Di/2pi/2 + ^) 2^Z}i/2 

Next, we use ffTTj) to express ^(x°) in terms of the initial data (11881) . This leads to 

y\x^) = M-^D-^'^^ X A(x°), y{x^) = SoA(x°). (190) 

As seen from (11851) . (I186p . and (I190p the action of the position and spin operators on a Proca field A 
mixes the components of the field. It is clearly more complicated than the action of the corresponding 
operators for the KG fields [H]. 

Next, we evaluate the action of the momentum, angular momentum, and helicity operators on 
A. Because Pq and p commute, in view of (I148p and (1591) . we have po = U'^'P'^U^.o. This in turn 
implies [po A](x°) = p A(x°) for all x° G M. Furthermore, using (11850 . (I186p . and (11900 . we can show 
that the angular momentum operator (in units of K) L := /i~^Xo x po acts on A according to 

^o(x°) = h-\^ X p)AO(x°) - y\x% ^ ^ .^gn 

^(a;0) = {;i-ixxp + S}A(x")-y(xO), ' ^ ' 

where ^{x^) := [Lyl](x°), and y{x^) := [soA](x°) is given in (I190p . Therefore, unlike the position 
Xq and spin Sq operators, the (linear) momentum po and the total angular momentum M := L + Sq 
operators have the same expressions as in nonrelativistic QM0 Similarly, we can compute the action 
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Note that the corresponding particle in nonrelativistic QM is described by A i.e., we should set = 0. 
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of the helicity operator on A. In view of fllQOp . we find: [(po ■ Sq) = (O, (p ■ S)A(x°)), which 

is in complete accordance with our previous result fllOOl) . 

Again, to the best of our knowledge, the actions of the position, spin, and other observables on a 
(covariant) Proca field A have not been previously given. Earlier works on the subject p!9l [20l [2T| [22] 
calculate the position and spin operators that act either on the six-component fields, i.e., in the 
Hilbert space /C, or on the Bargmann-Wigner's second-rank 4-spinors. One can use the unitary 
transformation p : IC ^ H' to compute the position (Xq := p^^Xgp) and spin (So := p^^S'p) 
operators acting in /C. The result is 

- " " 2A.I)V2(z)i/2 + ^) ^ ^) + " ^ " MD{Dy^ + M)^ /' ^'^'^ 

So = S f,/, ~^ ^(^ ■ S) + ^ ai {(j? • S)(i^ X S) + (j^ X S)(j^ • S)} . (193) 

These are exactly the position and spin operators that are obtained by Case in [T6] . 

Following the treatment of the KG fields in [271 128] we can identify the coherent states of Proca 
fields with the eigenstates of the annihilation operator a := (xo + ^fc'^Po), where G M. 

Because both xq and po commute with the chirality (ss) and spin (S12) operators, so does a. Hence, 
we can introduce a set of coherent states with definite chirality and spin. The corresponding state 
vectors |z, e, s) are defined as the common eigenvectors of a, S3 and S12, i.e., a|z,e, s) = z|z, e, s) 
and S3|z,e, s) = e|z,e, s), Si2|z,e, s) = s|z,e, s), where z G C^, e G {— ,+} and s G {—1,0,+!}. We 
can studied these coherent states and found that they have essentially the same properties as the 
coherent states of the KG fields that we explored in [28] . 



8 Probability Density for Spatial Localization of a Field 

We may employ the procedure outlined in [6] to find the probability density for the spatial localization 
of a Proca field. As in nonrelativistic QM, we identify the probability of the localization of a Proca 
field y4 in a region V C M^, at time to = x^/c, with 

Py=! ci^x lin^Af , (194) 
Jv 

where IIx is the projection operator onto the eigenspace of xo with eigenvalue x, i.e., 

nx = E E I4^'^^)(^^"^^l' (195) 

e=± s=0,±l 

II ■ II ^ := ((■,■)) is the square of the norm of 7i, and we assume ||y4|| = 1. Substituting fll95p in (11941) 
and making use of (I150p and (I15ip . we have Py = Ylit=±Ylis=o ±1 |/(e)'S,x)p. Therefore, in 

light of (11680 . fl52]) . and (1551) . the probability density is given by 

Q{xl^):=Y^ |/(e,s,x)r = -^{|HA(x°,x)p + |ilD-i/2A(x°,x)|n. (196) 

e=± s=0,±l 

For a position measurement to be made at time t = x^/c, we have the probability density 

^(x°,x) = ^ ||ilA(x°,x)p + |ilD-i/2A(x°,x)p} = ^ {|ilA(x)p + \iiA,{x)\'} . (197) 
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We can use the method discussed in Section |5]to introduce a current density such that = g 
(See also [50].) This probabihty current density turns out to have the form 



where 3? means "the real part of, T''(a;) = and, for all i G {1, 2, 3}, 



198) 



= ^5R<!(ilA(x))^ 



V • ilA(a;) 



M 



[^A(x)" ■ V] ( ^^iiAix) ) + 



(ilA,(x))^ 



V ■ ilAc(x) 



DV2 



M 



- [ilA.(x)" ■ V] ( ^j^iiA^ix] 



(199) 



One can easily show that J^{x) is neither a four- vector nor a conserved current density. 

Although the above discussion is based on a particular choice for the parameters {a}, namely 
{a} = {1}, it is generally valid. This is because the position wave functions / and the corresponding 
probability densities do not depend on the parameters {a}. Therefore, if we are to compute the 
probability density q^^-^ of the spatial localization of a Proca field A G ^^{a} with the position operator 
being identified with xo{a} for {o} ^ {!}, we have, for a measurement made at to = ^^o/c. 



Q^^^{xU) = ^ {|UA'^„^(a:°,x)p + lilD-^/^A'^^- 



}l4x)p 



(200) 



where A'|^| := ^ i^^A and U^a} '■ 'H ^ 7i{a} is given by fll53p . We can compute A'^a} using fl65 
(11201) . and (11531) and use the result to obtain 



{|il+,+A(x)|2 + |il_,+A(a;)p + |il+,+A,(x)|2 + |H_,+A,(x) 
+23? [(il+,+A(x))* ■ il+,+Ae(x) - (il_,+A(x))* ■ il_,+Ae(x)]} . 



Am 



(201) 



Again we can use the method of Section [5] to compute a probability current density J7|^| such that 
J'^^y = g^^y This leads to a complicated expression that we do not include here. Similarly to J'^, 
we expect J7|^| to be neither covariant nor conserved. 

The non-conservation (respectively non-covariance) of the probability current density i7|^| raises 
the issue of the non-conservation (respectively frame-dependence) of the total probability: 



Jm.3 



(202) 



This would certainly be unacceptable. The situation is analogous to that of the KG fields. /!P{a} 
is indeed a frame-independent conserved quantity, thanks to the covariance and conservation of the 
current density J^^y and the identity 



£x ^^„j(a;°,x) = / d^x Jf„}(x°,x), 



(203) 



which follows from (11281) . (12011) and the fact that il^^e/ are self-adjoint operators acting in H. 
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Combining fl202D and fl203|) . we have 

^{,}= [ d=^x J{'„^(x°,x). (204) 

This relation imphes that although the probability density g^^^ is not the zero-component of a con- 
served four-vector current density, its integral over the whole space that yields the total probability 
(12021) is nevertheless conserved. Furthermore, this global conservation law stems from a local con- 
servation law, i.e., a continuity equation for a four- vector current density namely J'^^y 

9 Gauge Symmetry Associated with the Conservation of the 
Total Probability 

In this section we explore a global gauge symmetry that supports the conservation of the total 
probability or its local realization as the conservation of the current density J^a}- determine the 
nature of this symmetry, we recall that the conserved charge associated with any conserved current is 
the generator of the corresponding infinitesimal gauge transformations [51]. We use the Hamiltonian 
formulation to obtain these transformations. The procedure we follow mimics the one presented in 
[6] for the KG fields. 

The Lagrangian L for a Proca field A and the corresponding canonical rnomenta n(x), n(x) 
associated with A(x) := A{x^,x.) and A*{x^,x.) are respectively given by [3]: 
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L := - j^^ |-F^,(x)*F^-(x) + ^M^(x)*A'^(x)| , (205) 

n°(x) := = 0, (x) := 4^ = -F°^(x)* = -E^(x)*, (206) 

Mo(x) ()A(x) 

n^x) := = 0, ff (x) := -i^ = -F°Xx) = -E\^). (207) 

Mo(x)* oAi{x.)* 



The fact that 11'' anc 



primary constraints: 



n'' vanish show that the Proca system is a constrained system. There are two 



E 

$1 := n°(x) ^ 0, $2 := n°(a;) ^ 0. (208) 

Solving for the velocities Ai and A* in (12061) and (120 7p . and using the Hamiltonian 

Ho:= I rf'x{n,(x)Il^(x)-ff(x)9,Ao(x)-ff(x)9,Ao(x)* + ^%(x)A°(x)*+ 

MMi(x)AXx)* + 9,A,(x)9M^(x)* - d,A.j{^)& A\^y] , (209) 



we obtain the so-called total Hamiltonian [52j: Ht '■= Hq + u^^j, where ,j = 1, 2 are two unknown 
coefficients. Using the Poisson bracket 



_M^(x) (5HA'(x) M^(x)5H'^(x) M^(x)*5HA'(x) M^(x)* 5Ha'(x) 



(210) 



^''Throughout this section we suppress the x^-dependence of the fields for simphcity. 

^"'^ Following Dirac's notation [52] . we write the constraints as weak equations with the weak equality symbol 
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of the observables T and we can easily show that the dynamical consistency of the primary 
constraints fl208p . i.e., $j = {$-,-, ifr}p ~ 0, results in the following secondary constraints, 

$3 := VAo(x) + (x) ^ 0, $4 := >f^^o(x)* + 9iff(x) ^ 0, (211) 

and that there is no other secondary constraint. 

It is not difficult to show that the matrix ^/(y, z) := {$j(y), $j/(z)}p is nonsingular. Hence we 
have a theory with four second-class constraints and can apply Dirac's canonical quantization that 
uses the Dirac bracket |52j : 

{^,^}^ := / d'y [ rf3z{^,<|.,(y)}p^^-^'(y,z){<|.,,(z),^}p, (212) 

where (y, z) is the inverse of ^^'(y, z). Computing the latter and its inverse, we find^^'^' (y,z) : = 
^n's^^y-z), and 

{J',g}^ = {J',g}p- [ d'y{J^,^,{y)}p^"'{^y{y),g}p, (213) 

where = '#^4 = _^3i ^ _^42 ^ ^-2^ ^^^^ c^jf ^ g ^^j^^^ ^-^^ ^^^^^ y,g^ Further details of the 
constraint quantization of Proca system can be found in [531 EH ES] . 

In terms of the canonical phase space variables (^,11) and [A*,]!), the total probability (12040 
takes the form 

= ^ I (n^x) \q+,oD~'/'a{^)] - M^y [e+,o/^-^/^n(x) 

+t (ff (x) [e+,oA(x)]^ - Ai^y [e+,on(x)]^)} , (214) 

where we have made use of fll28p . (12061) . and (I207p . Now, we can obtain the infinitesimal symmetry 
transformation, 

A^A + 6A, (215) 

generated by ^{a} using 

M(x) = {A(x),^{„l}^50, (216) 

where 5(f) is an infinitesimal real parameter. In view of (I206D . (I207l) . (I2T0|) . (|2T3|1 - (l2T6|l . (|29l). (|33|). 
and (|82|) . we have 

M°(x) = {L° C + L° } A°(x), (217) 

(5A(x) = -ise {e+,oC + e_,o} a(x), (218) 



where 69 := K6(j)/{2M) and 6±^o are given by (1^ . We may employ (11001) to express (I217p and (j^ 



as 

M(x) = -^5^^ {79+,oC + 79_,o} A(x), (219) 
where 'd^^ : 'H{a} — ^ "^{a} is defined by 

^,,0 := [Lt - Ll]S)' + L:S) + Ll (220) 
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In view of (12191) . the symmetry transformations (12151) are generated by the operator 'i?+,oC + 
i9_^Q. We can easily exponentiate the latter to obtain the corresponding non-infinitesimal symmetry 
transformations, 

A ^ e-K^+.oC+i9..o)^^ (221) 
where ^ G M is arbitrary. In terms of the (e, /;,)-components A^^h of A, (12211) takes the form 

^ = E E ^^.'^ ^ E E e"""^''^^^,^' (222) 

e=± h=0,±l e=± h=0,±l 

where we have made use of = 1 and S)^ = S^. 

Similarly to its spin- zero counterpart [6], as seen from (I22ip and (I222p . the gauge group@ G^a} 
associated with these transformations is a one-dimensional connected Abelian Lie group. Therefore, 
it is isomorphic to either of U{1) or R"*", the latter being the noncompact multiplicative group of 
positive real numbers |56]. 

We can construct a faithful representation of the group G{a.} using the six-component represen- 
tation A = y4_|_ _i, yl+ 0) ^-,+1) ^-,-1) ^-,o)"^ where C and are, respectively, represented by 
S3 and S12. In this representation a typical element of G{a} takes the form 

ga{9) := diag (e''"+'+^^ e-'^+-''^, e-''^+'°^ e'"--+^^ e'°—'^, e'"-"^) . (223) 

This expression suggests that the gauge group G^a} is a subgroup of f/(l)®^, the latter being the 
direct product of six copies of U{1). It is not difficult to show that G^a} is a compact subgroup of 
this group and consequently isomorphic to U{1) if and only if all the parameters a^^h are rational 
numbers, otherwise G^a} is isomorphic to M^. 

Clearly, the G{a} gauge symmetry associated with the conservation of the total probability is a 
global gauge symmetry. Similarly to its spin- zero counterpart [6] the local analog of this global gauge 
symmetry is different from the usual local Yang-Mills-type gauge symmetries. 

10 Conclusion and Discussion 

In this article we have used the methods of pseudo-Hermitian quantum mechanics to devise a com- 
plete formulation of the relativistic quantum mechanics of the Proca fields that does not involve 
restricting to the positive-energy solutions of the Proca equation. In particular, we have constructed 
the most general physically admissible inner product ((■,-)){a} on the solution space of the Proca 
equation. Up to a trivial scaling, this inner product involves five real parameters that we collectively 
denote by a. For all the values of a the inner product ((-, ■)){a} is positive-definite and relativisti- 
cally invariant. It also renders the generator of the time-translations, i.e., the Hamiltonian h, and 
the helicity operator self-adjoint. The quantum system associated with the Proca fields may be 
represented by a Hilbert space 'H{a} defined by the inner product ((■, ■)){a} and the Hamiltonian h. 

^^Here we identify the gauge group with its connected component that includes the identity and is obtained by 
exponentiating the generator i^+.qC + i9^^o- 
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Using the unitary equivalence of this covariant representation with that of the noncovariant Foldy- 
representation, which arises quite naturally in our formulation, we constructed relativistic position, 
momentum, angular momentum, spin and helicity operators acting in 7i{a}, and the localized states of 
the Proca field. Furthermore, we introduced the position wave functions and used them to construct 
a probability current density which turned out to be neither conserved nor covariant. We resolved 
the apparent inconsistency of this observation with the physical requirement of the conservation and 
frame-independence of the total probability using a conserved and covariant current density. The 
global conservation of the total probability is supported by the local conservation of this current den- 
sity. The latter is linked to a previously unnoticed global gauge symmetry of the Proca field with an 
Abelian gauge group. In Ref. [29], the authors have also attempted to give a spin-one generalization 
of our treatment of the Klein-Gordon fields. However, because of various self-imposed restrictions 
they obtain a one-parameter subfamily of the inner products ((■, ■)){a}@ This is related to the fact 
that the authors of [29] use the special Foldy transformation fl65|) (see Eq. (35) of [29]) to find the 
metric operator acting in the Hilbert space Ti' . Furthermore, in trying to impose the condition of 
the Lorentz-invariance of the metric operator, they transform the operator f/+ to a metric operator 
acting in the Foldy representation and demand that the latter commutes with the generators of the 
Poincare group in this representation. This does not seem to be well-justified, because as seen from 
(HOll . the Hilbert space of the Foldy representation is just the direct sum of two copies of 7i. Hence, 
the metric operator associated with this representation is just the identity operator, not the one 
given by Eq. (39) of [29j. As seen from (11171) . this is the Foldy transformation which depends on 
the choice of the unknown parameters appearing in the operator fjj^, not the metric operator of the 
Foldy representation. We would also like to stress that the analysis of [29] does not include the 
construction of the observables of the system or any treatment of its physical aspects such as the 
notorious problem of the probabilistic interpretation of the quantum mechanics of Proca fields. 

Apart from the historical importance of the subject, the present work is mainly motivated by 
the close analogy of the Proca and Maxwell fields. Performing the zero-mass limit of our results 
in an appropriate manner should lead to a consistent quantum mechanical treatment of individual 
photons. As it is to be expected, there are subtleties in performing this limit. Nevertheless, we have 
been able to make some progress toward solving the problem of the construction of the Hilbert space 
and observables for the photon. We plan to report the results in a separate article. 
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